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Search for Event Rate Modulation in XENONIOO Electronic Recoil Data 
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We have searched for periodic variations of the electronic recoil event rate in the (2 — 6) keV energy 
range recorded between February 2011 and March 2012 with the XENONIOO detector, adding up 
to 224.6 live days in total. Following a detailed study to establish the stability of the detector and 
its background contributions during this run, we performed an un-binned profile likelihood analysis 
to identify any periodicity up to 500 days. We find a global significance of less than 1 a for all 
periods suggesting no statistically significant modulation in the data. While the local significance 
for an annual modulation is 2.8 a, the analysis of a multiple-scatter control sample and the phase 
of the modulation disfavor a dark matter interpretation. The DAMA/LIBRA annual modulation 
interpreted as a dark matter signature with axial-vector coupling of WIMPs to electrons is excluded 
at 4.8 a. 

PACS numbers: 95.35.+d, 14.80.Ly, 29.40.-n, 
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The XENONIOO experiment [1] is designed to search 
for dark matter in the form of Weakly Interacting Mas¬ 
sive Particles (WIMPs) [2] by detecting WIMP-induced 
nuclear recoils (NRs) with a liquid xenon (LXe) time pro¬ 
jection chamber. The resulting event rate in any dark 
matter detector is expected to be annually modulated 


due to the relative motion between the Earth and the 
dark matter halo of the Milky Way [3]. The modula¬ 
tion of the low energy (low-E), (2 — 6)keV, event rate 
in the DAMA/LIBRA experiment [4] is currently the 
only long-standing claim for a positive dark matter de¬ 
tection. Under typical astrophysical and particle physics 
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assumptions, this claim is however challenged by the non¬ 
observation of WIMP-induced NRs of several other ex¬ 
periments using different target materials and detector 
technologies [e.g. 5-7], most with considerably lower ra¬ 
dioactive backgrounds. 

An alternative explanation is that the DAMA/LIBRA 
annual modulation is due to electronic recoils (ERs) 
from WIMPs which have axial-vector couplings to elec¬ 
trons [9, 10]. The stable performance of XENONIOO over 
a period of more than one year offers the opportunity to 
test this hypothesis with a different detector operated for 
the first time in the same underground site, the Labora- 
tori Nazionali del Gran Sasso (LNGS), Italy. 

For this analysis we use the 224.6 live days of 
XENONIOO dark matter data accumulated from Febru¬ 
ary 28, 2011 to March 31, 2012, previously used to search 
for spin-independent [5] and spin-dependent [11] WIMP- 
induced NRs as well as for axion-induced ERs [12] and 
a comparison with DAMA/LIBRA using the average ER 
rate [10]. 

The ER energy and uncertainty therein is inferred from 
the prompt scintillation light signal (SI), as in [12], us¬ 
ing the NEST model (v0.98) [13] fit to independent light 
yield calibration measurements [14, 15]. The overall un¬ 
certainty on the ER energy scale is dominated by the 
spread of the low energy measurements in [14, 15] and is 
estimated to be 14% at 2 keV and 9% at 6 keV. 

We use the same SI range of (3 — 30) photoelectrons 
(PE) as in [5, 16], but divided into two ranges. The low- 
fi range (3 — 14) PE corresponds to (2.0 — 5.8) keV and 
thus covers the energy interval where the DAMA/LIBRA 
experiment observes a modulation signal. The higher 
energy range, (14—30) PE, corresponds to (5.8—10.4) keV 
and is used as a sideband control sample. 

Low-E single-scatter events in the 34 kg fiducial mass, 
as expected from dark matter interactions, are selected 
using the same cuts as in [5]. While these cuts were 
defined to select valid NR events, they also have high ef¬ 
ficiency for ERs [12], and result in 153 events distributed 
in time as shown in Fig. 1(f). The cut acceptances in 
the energy ranges considered here have been derived fol¬ 
lowing the procedure in [16] using ER calibration data 
( 60 Co and 232 Th) taken on a weekly basis. The time 
variation of the acceptance, shown in Fig. 1 (e), is incor¬ 
porated in the analysis by linearly interpolating between 
the data points. We have verified that our conclusions re¬ 
main unaffected when adopting different methods of cut 
acceptance interpolation in time. 

The design of XENONIOO incorporates many sensors 
of various types to monitor the long-term stability of the 
detector. A total of 15 parameters were investigated, of 
which a subset with the highest potential impact on de¬ 
tector signals is shown in Fig. 1 (a-d). The absolute pres¬ 
sure of the gas above the LXe has a mean value of 2.23 bar 
with a maximum variation of 0.02 bar over the entire pe¬ 
riod (Fig. 1 (a)). The temperature sensors located at var- 
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FIG. 1: Temporal evolution of the relevant XENONIOO de¬ 
tector parameters studied for this analysis. The dashed blue 
lines indicate a detector maintenance period, (a-c) Xe pres¬ 
sure, LXe temperature and LXe level, (d) Radon level in the 
34 kg LXe fiducial mass, as measured via in situ alpha spec¬ 
troscopy. (e) Average cut acceptance in the low-E range of 
(2.0 — 5.8) keV, as derived from weekly ER calibrations, (f) 
ER event rate in the 34 kg fiducial mass for single-scatters in 
the low-E range. 
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ious positions within the detector exhibit a mean value 
varying from —86.5°C (in the Xe gas) to — 91.6°C (bot¬ 
tom of the LXe) with a maximum variation of less than 
0.17 °C for each sensor. The ambient temperature in the 
XENONIOO room has a mean value of 20.7°C with a 
maximum variation of 3.7°C (Fig. 1 (b)). The LXe level, 
monitored by two capacitive sensors, shows a maximum 
variation of 0.22 mm during the entire period (Fig. 1 (c)). 

To identify potential correlations between detector pa¬ 
rameters and ER rate, we calculate the linear (Pearson) 
and non-linear (Spearman-Rank) correlation coefficients 
for the two energy ranges studied, and for both single¬ 
scatter and multiple-scatter events. The latter are de¬ 
fined as events with a single-scatter in the fiducial re¬ 
gion plus an additional SI coincident signal in the LXe 
veto. The 99 kg LXe veto has an energy threshold of 
~100keV, thus multiple-scatter events are dominated by 
high-energy scatters from 7 rays [1, 19]. Of all the param¬ 
eters studied, two were found to give a non-correlation 
p-value smaller than 0.001. The first parameter is the 
LXe level, which shows a negative linear and non-linear 
correlation with the low-E single-scatter rate. The second 
parameter is the Xe gas temperature, which shows a neg¬ 
ative linear correlation with the low-E multiple-scatter 
rate. As expected, the LXe level and gas temperature 
were also found to be correlated with each other and with 
the room temperature. A change in the LXe level, gas 
pressure and temperature can potentially affect the ob¬ 
served size and width of the secondary scintillation signal, 
S2, which is a measure of the ionization electrons liber¬ 
ated in the interaction. The overall observed variation of 
the S2 signal is less than 5% [17], while the majority of 
events have S2 > 1000 PE, much larger than the trigger 
threshold of 150 PE. Consequently, a detailed inspection 
of the S2-dependent cuts shows that their performance 
is unaffected. Hence the correlation with event rate is 
possibly a coincidence and, regardless, does not impact 
our statistical analysis for periodicity described below. 

The impact of decaying radioactive isotopes on the low- 
fi ER rate is also considered for this analysis. These 
sources can be subdivided into external sources of 7 - 
radiation from peripheral materials and /3-radiation from 
the decay of radioactive Rn and 85 Kr distributed in the 
LXe volume. 

Of the relevant external 7 -sources in the detector and 
shield materials, only 60 Co (ti / 2 = 5.27y) decays on a 
timescale sufficiently short to potentially cause an ob¬ 
servable change in the event rate during the time period 
of this study. However, the decrease in activity is found 
to reduce the single-scatter low-E ER rate by less than 
1% of its average value, based on a Monte Carlo (MC) 
simulation using the measured activity level from [19]. 
Hence we assume the external 7 -background to be con¬ 
stant for this analysis. 

The short-lived isotopes 222 Rn and 220 Rn are con¬ 
stantly produced as part of the primordial 238 U/ 232 Th 


decay chains and are present in the air of the room and 
shield cavity, as well as inside the LXe due to emanation 
from inner surfaces. Radon decays outside the detector, 
measured by commercial Rn monitors in the room, con¬ 
tribute negligibly to the event rate in the fiducial mass 
since the emitted radiation is absorbed by the shield and 
outer detector materials. The concentration of Rn and 
subsequent decay products dispersed in the LXe is con¬ 
tinuously monitored via examination of both a-decays 
and /3-q delayed coincidence events [18]. This analy¬ 
sis shows that 222 Rn from the 238 U chain is uniformly 
distributed in the volume while 220 Rn from the 232 Th 
chain is negligible. The time-variation of the 222 Rn level 
is shown in Fig. 1 (d) and exhibits a specific activity of 
(63 ± 1) /iBq/kg. This level corresponds to a low-E ER 
contribution of (1.11 ±0.02) events/(keV • tonne • day) as 
determined by MC simulation [19]. The 8.5% fluctuation 
of the 222 Rn level corresponds to a less than 2% variation 
of the average rate and is thus negligible compared to the 
observed rate fluctuation of 51% shown in Fig. 1 (f). In 
addition, no time correlation is found by calculating the 
linear and non-linear correlation coefficients between the 
low-E ER rate and the Rn level. Therefore the evolution 
of the 222 Rn level in time is not included in the statistical 
analysis below. 

The other internal contamination, 85 Kr, is also present 
in air. The concentration of nat Kr in the LXe during 
the period studied here was determined on November 
17, 2011 to be (14 ± 2) parts per trillion using the rare 
gas mass spectrometer (RGMS) method [5, 20]. How¬ 
ever, it became evident after the end of the run that a 
small air leak in the Xe gas purification system had al¬ 
lowed Rn and Kr atoms to diffuse into the LXe. The 
leakage rate into the sensitive volume was estimated 
from a study of the time correlation between the ex¬ 
ternal and internal concentrations of 222 Rn [18], includ¬ 
ing three RGMS measurements of nat Kr spread over the 
course of several months during the following run. As¬ 
suming a constant nat Kr concentration in air, the lin¬ 
ear increase in time of nat Kr in the LXe was found 
to be proportional to the integrated number of addi¬ 
tional 222 Rn decays due to the air leak. The linear 
increase of the single-scatter ER rate from 85 Kr has a 
slope K = (2.54 ± 0.53) x 10 -3 events/(keV • tonne • 
day)/day assuming a 85 Kr/ nat Kr ratio of 2 x 10 -11 [20]. 
This time-dependent background results in an expected 
total increase of (0.10 ± 0.02) events/day at low-E over 
the course of one year, which is taken into account in the 
following statistical analysis. 

To determine the statistical significance of a periodic 
time dependence in the event rate, we implement an 
un-binned profile likelihood (PL) method [21], which 
incorporates knowledge of the time variation of detec¬ 
tor parameters and radioactive backgrounds as described 
above. The event rate for a given energy range is de- 
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Frequency [1/Days] 

0.1 0.02 0.01 0.002 0.001 



FIG. 2: The expected mean (solid lines) and central 

68.3% region (shaded bands) of — 2log(£o/£i) as a func¬ 
tion of period for simulated data with a fixed average rate 
C = 6.0 events/(keV • tonne • day), linear increase in rate 
K = (2.54±0.53) x 10~‘ ! events/(keV ■ tonne ■ day)/day, am¬ 
plitude A = 2.7 events/(keV • tonne • day), and three periods 
P [days]. Uncertainties on all parameters are taken into ac¬ 
count. The horizontal local significance lines are derived from 
the null hypothesis tests described in the text and shown here 
for comparison to Fig. 3. 

scribed by 

f{t) = e(t) (c + Kt + Acos ^2-7 t — > (!) 

where e is the corresponding average cut acceptance, 
interpolated from the measurements described above, 
C is the constant component of the event rate, Kt is 
the linearly increasing contribution from 85 Kr, and A is 
the modulation amplitude with period P and phase <j>. 
Eq. (1) is then normalized to take into account the time 
distribution of the dark matter data used for the analysis 
here, and thus becomes the probability density, f(t), of 
observing an event occurring at time t, in days relative 
to January 1, 2011. The null hypothesis, no periodicity, 
is given by Eq. (1) with A = 0. 

The likelihood function used in the PL method is 

C = /(£) j Poiss (nliVexp^)) C e C K C E , (2) 

where n and N exp (E) are the total number of observed 
and expected events and E is the energy in keV. Nuisance 
parameters corresponding to the uncertainties in e, K, 
and E are constrained by the Gaussian penalty terms, 

. Ck , and C E , respectively. These penalty terms have 
widths a e defined by the statistical errors of the accep¬ 
tance as determined by weekly calibration measurements, 
<j k = 0.53 x 10~ 3 events/(keV • tonne • day)/day, and 
a E taken from Fig. 2 of [12], respectively. The maximum 
profiled likelihoods are denoted by £o(Co) for the null 
hypothesis and by £i(Ci,A,<j>) for the periodic hypoth¬ 
esis. 

The significance of a particular period, for example 
one year, is referred to as the local significance. The 
corresponding test statistic is the log-likelihood ratio, 
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FIG. 3: — 21og(£o/£i) as a function of modulation period 

for single-scatters (SS) in the low-E region (top), multiple- 
scatters (MS) in the low-E region (middle) and single-scatters 
(SS) in the higher energy region (bottom). The phase is un¬ 
constrained. 

— 21 og(£o/£i), which quantifies the incompatibility be¬ 
tween the null and periodic hypotheses. MC simulations 
show that this test statistic is well-described by a y 2 - 
distribution with two degrees of freedom. When search¬ 
ing for a modulation signal across a range of periods, the 
global significance, that is the maximum of the local test 
statistics in the range, should be referenced. The local 
and global significances quoted are both one-sided. 

Simulated data were used to assess the discovery po¬ 
tential of the PL analysis to periodic components in the 
single-scatter data at low-E. Several sets of 153 simu¬ 
lated events were generated by drawing from the same 
live-time distribution as the actual data while varying 
the nuisance parameters according to their constraints in 
Eq. 2, and assuming the periodic hypothesis with a fixed 
period, amplitude and average rate. The expected signif¬ 
icance is shown in Fig. 2 for three periods with an ampli¬ 
tude of 2.7 events/(keV • tonne • day) and average rate of 
6.0 events/(keV • tonne • day), selected to facilitate com¬ 
parison with the best-fit results below. The minimum 
period considered is 7 days, since the cut acceptance is 
derived from weekly calibration measurements. The res¬ 
olution on the reconstructed period becomes worse with 
increasing period, evident from the broadening of the 
peaks and a characteristic plateau for periods > 500 days. 
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Hence the study of the data in Fig. 3 was limited to peri¬ 
ods between 7 and 500 days. Adding the previous 100.9 
live days of data [22] to this analysis does not consid¬ 
erably increase the significance of the study due to its 
higher background rate from 85 Kr and the uncertainty 
therein. 

In addition to the un-binned PL analysis, a x 2_ test 
following [23] and a Lomb-Scargle (LS) periodogram [24] 
were carried out using binned data. For both tests, a 
strong binning dependence of the result is observed. This 
dependence, as well as the unavoidable information loss 
when using any bin-dependent method, limits the power 
of these tests compared to the un-binned PL analysis. 
This fact must be taken into account when using the data 
in Fig. 1 (f) for further analysis. Nevertheless, the local 
and global significances are in agreement with the results 
of the PL analysis and the tests provide a consistency 
check. 

WIMP interactions in the LXe are expected to produce 
single-scatter events. The PL spectrum of the single¬ 
scatter data covering the DAM A/LIBRA energy region 
(2.0—5.8 keV) is shown in the top panel of Fig. 3. A rise in 
significance is observed at long periods with a local signif¬ 
icance of 2.8 a at one year and a global significance below 
1 a for all periods. MC simulations with P = 100 days in 
Fig. 2 show that the rise of significance at large periods 
in the measured data is not an artifact of the statistical 
method. 

Low-E multiple-scatter events are used as a 
background-only control sample. The PL spectrum 
(middle panel of Fig. 3) shows a rise in significance at 
long periods, similar to that for single-scatters, with 
a local significance of 2.5 a at one year and a global 
significance below 1 a at all periods. 

As WIMPs are expected to produce signals primar¬ 
ily at low-E, the higher energy range (5.8 — 10.4 keV) 
is used as a sideband control sample. In addition, 
DAM A/LIBRA did not observe a modulation above 
6 keV. The PL spectrum (bottom panel of Fig. 3) shows 
no prominent rise in significance at long periods, in con¬ 
trast to that seen at low-E, and the local significance is 
1.4 a at one year. 

The best-fit parameters and uncertainties are deter¬ 
mined from PL scans. For an assumed annual modula¬ 
tion signal (fixing P = 365.25 days) in the low-E single 
scatter data, we obtain C\ = (5.5 ± 0.6) events/(keV • 
tonne • day) (for reference, Cq = 6.0 events/(keV • tonne 
• day)), A = (2.7 ± 0.8) events/(keV • tonne • day), and 
</> = (112 ± 15) days, peaked at April 22. Fig. 4 shows 
the corresponding confidence level contours as a function 
of modulation amplitude and phase. The simulations 
in Fig. 2 show that the rise in significance at long peri¬ 
ods in the low-E single- and multiple-scatter data could 
be explained by a modulating component with a period 
>300 days. However, the best-fit phase disagrees with 
the expected phase from a standard dark matter halo 



Phase [Days] -21og(£,/£ m J 


FIG. 4: The XENON100 best-fit, 95% and 99.73% confi¬ 
dence level contours as a function of amplitude and phase 
relative to January 1, 2011 for period P = 1 year. The 

expected DAMA/LIBRA signal with statistical uncertainties 
only and the phase expected from a standard dark matter 
(DM) halo are overlaid for comparison. Top and side panels 
show — 2 \og(C\/ Cmax) as a function of phase and amplitude, 
respectively, along with two-sided significance levels. 


(152 days) at a level of 2.5 er based on the ID PL scan 
as shown in top panel of Fig. 4. Furthermore, the rise in 
significance at long periods is evident in both single- and 
multiple-scatter data, also disfavoring a WIMP interpre¬ 
tation. Allowing the parameter K to float freely to un¬ 
physical negative values, given the measured 85 Kr level, 
decreases the significance of large periods and strength¬ 
ens the exclusion limit discussed below. 

The XENON100 data can constrain the dark mat¬ 
ter interpretation of the annual modulation observed 
by DAMA/LIBRA, as shown in Fig. 4, for certain 
models producing ERs. Such constraints were pre¬ 
viously imposed using the average ER event rate in 
XENON100 [10]. Here we use the full time-dependent 
rate information to directly compare with the expected 
DAMA/LIBRA annual modulation signal in our detec¬ 
tor. The expected SI spectrum in XENON100 is de¬ 
rived from the DAMA/LIBRA residual modulation spec¬ 
trum (Fig. 8 in [4]) following the approach described 
in [10], assuming the signals are from WIMP-electron 
scattering through axial-vector coupling [9, 10]. The ex¬ 
pected annual modulation amplitude in the low-E range 
in XENON100 is then calculated as (11.5 ± 1.2(stat) ± 
0.7(syst)) events/(keV • tonne • day), with statistical 
uncertainty from the reported DAMA/LIBRA spectrum 
and systematic uncertainty from the energy conversion in 
XENON100. To compare this expected signal with our 
data, the phase (/ in Eq. (1) is set to (144 ± 7) days [4], 
constrained by an additional Gaussian term, in Eq. 2. 
The resulting PL analysis of our data disfavors the ex¬ 
pected DAMA/LIBRA annual modulation at 4.8 er. 

In summary, XENON100 has demonstrated for the 
first time that LXe dual-phase time projection cliam- 
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bers can be operated with sufficient long-term stability 
to enable searches for periodic signals for periods up to 
and greater than one year. The detector parameters in¬ 
vestigated were found to be very stable, and most show 
no correlation with the measured low-E (2.0 — 5.8 keV) 
single-scatter ER event rate. Although the LXe level and 
Xe gas temperature show a correlation with this rate, no 
evidence was found of a direct impact on the cut perfor¬ 
mance. A time varying cut acceptance and background 
from 85 Kr are included in the search for event rate mod¬ 
ulation. In the 224.6 live days of XENONIOO data taken 
over the course of more than one year, a rising signifi¬ 
cance at long periods is observed for low-E single- and 
multiple-scatter events with the most likely period be- 
ing >450 days. An explicit search for annual modula¬ 
tion in the ER rate gives a 2.8 a focal significance with 
a maximum rate on April 22 ± 15 days. This phase dis¬ 
favors an annual modulation interpretation due to the 
standard dark matter halo at 2.5 a. Furthermore, our 
results disfavor the interpretation of the DAMA/LIBRA 
annual modulation signal as being due to WIMP-electron 
scattering through axial-vector coupling at 4.8 er. 
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